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Abstract. The modified Newtonian dynamics (MOND) has been proposed as an 
alternative to the dark matter paradigm; the philosophy behind is that there is no dark 
matter and we witness a violation of the Newtonian law of dynamics. In this article, 
we interpret differently the phenomenology sustaining MOND, as resulting from an 
effect of "gravitational polarization" , of some cosmic fluid made of dipole moments, 
aligned in the gravitational field, and representing a new form of dark matter. We 
invoke an internal force, of non-gravitational origin, in order to hold together the 
microscopic constituents of the dipole. The dipolar particles are weakly influenced by 
the distribution of ordinary matter; they are accelerated not by the gravitational field, 
but by its gradient, or tidal gravitational field. 
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1. Introduction 

The observed discrepancy between the dynamical mass and the luminous mass of 
bounded astrophysical systems is generally attributed to the existence of an invisible 
form of matter, coined as the missing mass or dark matter. The nature of the dark 
matter particles is unknown but extensions of the standard model of particle physics 
provide a number of candidates pQ. The dark matter triggers the formation of large- 
scale structures by gravitational collapse and predicts the scale-dependence of density 
fluctuations. Simulations suggest some universal dark matter density profile around 
ordinary matter distributions [2j. An important characteristic of dark matter, required 
by the necessity of clustering matter on small scales, is that it should be cold or non- 
relativistic at the epoch of the galaxy formation. However, the dark matter hypothesis 
has some difficulties [3J at explaining naturally the flat rotation curves of galaxies, one 
of the most persuasive evidence for the existence of dark matter, and the Tully-Fisher 
empirical relation |4| between the observed luminosity and the asymptotic rotation 
velocity of spiral galaxies. 

On the other hand, the modified Newtonian dynamics (MOND) has been proposed 
by Milgrom [5J, [6], E] as an alternative to the dark matter paradigm. It imputes the 
mass discrepancy not to the presence of some additional non-baryonic matter, but to 
failure of the Newtonian law of gravity in an appropriate regime — a drastic change of 
paradigm. MOND involves a single parameter a being a constant acceleration scale, 
which delineates the specific MOND regime, corresponding to accelerations much smaller 
than a , from the Newtonian regime, for which the accelerations are much larger. Several 
relativistic extensions of MOND, assuming the existence of extra fields associated with 
gravity, besides the spin-2 metric field of general relativity, have been proposed [T3l fTi]. 
Such extensions have culminated in the scalar-vector-tensor theory of Bekenstein and 
Sanders [ElIIElin]. 

MOND has been very successful at fitting the flat rotation curves of galaxies, and 
at recovering naturally the Tully-Fisher relation (see [HJ [TOj, [TT], [12] for reviews). 
Intriguingly, the numerical value of a that fits the data is close to the Hubble scale, 
ao ~ cH Q . MOND may have also some observational problems: There are some 
counter-examples of galaxies where MOND does not seem to account for the observed 
kinematics [18J, and, most importantly, the mass discrepancy at the scale of clusters of 
galaxies is not entirely explained by MOND [T9l 120] . 

In this article and the next one [21] (hereafter paper II), we take the view that 
MOND does not represent a violation of the fundamental law of gravity, but, rather, 
provides us with an important hint on the (probably unorthodox) nature of the elusive 
dark matter. More precisely, we interpret the phenomenology behind MOND as resulting 
from an effect of gravitational polarization, of some cosmic fluid made of dipole moments, 
and representing a new form of dark matter. The dipole moments get aligned in the 
gravitational field produced by ordinary masses, thereby enhancing the magnitude of 
the field and yielding MOND. Such effect is the gravitational analogue of the usual 
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electrostatic effect, of polarization of a dielectic medium in an applied electric field |22j. 

In the present paper we imagine, as a model for the dipole, a doublet of particles, 
one having a positive gravitational mass, the other having a negative and opposite 
gravitational mass, and with both particles being endowed with positive inertial masses. 
The gravitational behavior involving masses of this type is governed by a negative 
Coulomb law — like masses attract and unlike masses repel [23J . As a result the dipole 
moment cannot be stable: Even if we neglect the repulsive gravitational force between 
the two particles, they will accelerate apart from each other in an exterior gravitational 
field produced by ordinary matter. We shall therefore invoke an internal force, of non- 
gravitational origin, between the two particles constituting the dipole, to bound them 
in a gravitational field. The MOND acceleration scale a will appear to be related to 
the properties of this internal "microscopic" force at short distances. We find that the 
motion of the dipolar particles violates the equivalence principle, and is driven by the 
tidal gravitational field of ordinary matter, rather than the gravitational field itself. In 
this sense the dark matter is only weakly influenced by the distribution of ordinary 
matter. 

Summarizing, in our approach the dark matter is described by a "digravitational" 
medium, which is subject to polarization in a gravitational field, and is otherwise 
essentially static (an "ether"). An alternative interpretation of this dark matter is by 
the gravitational analogue of a plasma in electromagnetism, i.e. composed of positive 
and negative gravitational masses, and oscillating at the natural plasma frequency. In 
a gravitational field the mean position of the masses is displaced from equilibrium and 
the plasma acquires a dipolar polarization. The observational predictions of the present 
(non-relativistic) model are the same as for MOND. 

In paper II we shall propose a relativistic model of dipolar particles, based on an 
action principle in general relativity. This model will be consistent with the equivalence 
principle, and as a result the dynamics of dipolar particles, even in the non-relativistic 
limit, will be different from that of the present quasi-Newtonian model. The present 
paper and paper II provide two distinct models, both of them suggest a close connection 
between the phenomenology of MOND and some form of gravitationally polarized 
dipolar dark matter. 

Section [2] investigates the formal analogy between the MOND equation and the 
electrostatics of non-linear media. Sections [3] and 0] introduce a microscopic quasi- 
Newtonian description of the gravitational dipole moment. In Sec. |5] we derive the 
expression of the non-gravitational internal force in the MOND regime. In Sec. [6] we 
present an alternative though equivalent formulation of the dipolar medium in terms of 
a polarized gravitational plasma. 
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2. Analogy with electrostatics 

The MOND equation, in the variant derivable from a non-relativistic Lagrangian [13J, 
takes the form of the modified Poisson equation 



0; 



/'I — J 

a 



AwGp, (1) 



where p denotes the density of ordinary matter, g % = d^J is the gravitational field in the 
non-relativistic limit (so that a % = g % is the acceleration of ordinary matter), and U is the 
gravitational potential Jj| In Eq. (JTJ the Milgrom function p depends on the ratio g/a , 
where g — \g l \ is the norm of the gravitational field, and ao is the constant acceleration 
scale. The MOND regime corresponds to the limit of weak gravity, much below the 
scale cio, i.e. g <C a ; in this limiting regime we have p(g/a ) « g/a [HI EE E]- When 
g ^> a , the function p,(g/a ) asymptotes to one, and we recover the usual Newtonian 
law. Sometimes we shall consider the formal "Newtonian" limit g —>■ oo; however 
we always assume that the gravitational field is non-relativistic. Various forms of the 
function p have been proposed (c/. [211), but most of them appear to be rather ad hoc. 

Taking the MOND equation (CEJ) at face, we notice a striking analogy with the usual 
equation of electrostatics describing the electric field inside a dielectric medium, namely 
diD 1 = p e , where the electric induction D l is proportional to the electric field E l (at 
least for not too large electric fields): D % = p e eoE\ Here p e = 1 + Xe is the dielectric 
coefficient, and Xe denotes the electric susceptibility of the dielectric medium, which 
depends on the detailed microscopic properties of the medium (see e.g. [22] )■ In non- 
linear media the susceptibility is a function of the norm of the electric field, Xe{E) with 
E = \E t \. The electric polarization is proportional to the electric field, and is given by 
Ilg = Xe^oE 1 . The density of electric charge due to the polarization is p^° l = — diW e . 
Generically we have Xe > 0, which corresponds to screening of electric charges by the 
polarization charges, and reduction of the electric field inside the dielectric. 

In keeping with this analogy, let us interpret the MOND function p entering Eq. ([1]) 
as a "digravitational" coefficient, and write 

/^ = 1 + X, (2) 

where x would be a coefficient of " gravitational susceptibility" , parametrizing the 
relation between some " gravitational polarization" , say IT*, and the gravitational field: 

Since as we have seen the MOND function (T5]) depends on the magnitude of the 
gravitational field, p(g/a ), the same is true of the so defined gravitational susceptibility, 
x(g/ao), in close analogy with the electrostatics of non-linear media. Hence we expect 
that x should characterize the response of some non-linear digravitational medium to 

| Spatial indices i,j are raised and lowered using the Euclidean metric the summation convention 
is used throughout. 
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an applied gravitational field. The mass density associated with the polarization would 
then be given by the same formula as in electrostatics, 

Ppol = -a,ir . (4) 

With those notations Eq. ([1]) can be rewritten as 

Af/ = -47rG(p + p pol ). (5) 

In such rewriting of MOND, we see that the Newtonian law of gravity is not violated, 
but, rather, we are postulating the existence of a new form of matter, to be called dark 
matter, and which contributes in the normal way to the right-hand-side (RHS) of the 
Poisson equation (jSJ). The dark matter consists of polarization masses with volume 
density p pol given by (gj). 



3. Sign of the susceptibility coefficient 

For the moment we have restated the MOND equation in the form (j4])-(j5]) and proposed 
a formal interpretation. To check this interpretation, let us view the digravitational 
medium as consisting of individual dipole moments it 1 with number density n, so that 
the polarization vector reads 

rP = n7r\ (6) 

We suppose that the dipoles are made of a doublet of sub-particles, one with positive 
mass +m and one with negative mass —m, the masses which we are referring here 
are the gravitational masses of these particles, i.e. the gravitational analogue of the 
electric charges, m g = ±m. Clearly the exotic nature of this dark matter shows up here, 
when we suggest the notion of negative gravitational masses. If the two masses ±m are 
separated by the spatial vector d\ pointing in the direction of the positive mass, the 
dipole moment is 

7r* = m£f. (7) 

Let us further suppose that the two sub-particles are endowed with inertial masses which 
are positive, and given by m ; = m. The dipole moment thus consists of an ordinary 
particle, say (m i; m g ) = (m, m), associated with an exotic particle, [m- u m g ) = (m, —m). 

The ordinary particle will always be attracted by an external mass distribution 
made of ordinary matter; however, the other particle (m[,m g ) = (m, — m) will always 
be repelled by the same external mass. In addition the two sub-particles will repel each 
other. We see therefore that the gravitational dipole is unstable, and we shall invoke 
a non-gravitational force to sustain it. We also expect that the external gravitational 
field will exert a torque on the dipole moment in such a way that its orientation has 
the positive mass +m oriented in the direction of the external mass, and the negative 
mass — m oriented in the opposite direction. We thus find that 7r* and IT* should 
both point towards the external mass, i.e. be oriented in the same direction as the 
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external gravitational field g\ From Eq. fl3]) we therefore conclude that the susceptibility 
coefficient \i in the gravitational case, must be negative: 

X <0. (8) 

This corresponds to an " anti- screening" of the ordinary mass by the polarization masses, 
and enhancement of the gravitational field in the presence of the digravitational medium. 
This simply results from the fact that, in contrast to electrostatics, alike gravitational 
charges or masses always attract. The result (jSJ) is nicely compatible with the prediction 
of MOND; indeed we have < \l < 1 in a straightforward interpolation between the 
MOND and Newtonian regimes, hence — 1 < x < 0. The stronger gravitational field 
predicted by MOND may thus be naturally interpreted by a process of anti-screening 
by polarization masses. 



4. Microscopic model for the dipole 

To give more substance to the model, suppose that some interaction F % between the 
two constituents of the dipole is at work, and let the dipole moment be embedded 
into the gravitational field g % = diU . The equations of motion of the sub-particles 
(mi,m g ) = (m, m) and (m, — m), having positions x\ and x l _ respectively, are 

(fix 1 

m-^ = mg\x + )-F\x + -x^) 1 (9) 

m——= - mgHxJ) + F l (x + - x) . (10) 
at 2 

The internal force F l is proportional to the relative separation vector d l = x l + — x l _, 
namely 

F l = F-. (11) 

The norm of F % is a function of the separation distance, F = F(d) where d = \x\ — x l _\, 
and is expected to also depend on the magnitude of the gravitational field, g = \g l \. 
The force F l is assumed to be attractive, F > (hence it is non-gravitational). This 
force is indispensable if we are looking for configurations in which the constituents of 
the dipole remain at constant distance from each other. For simplicity, we suppose that 
the gravitational force between the two sub-particles, i.e. F g = —Gm 2 /d 2 , which is 
repulsive, is negligible or included into the definition of F. 

We introduce next the centre of inertial masses, x % = [x\ + x l _)/2, and transform 
the system of equations (!9l- (fT0l) into an equation of motion for the "dipolar particle" , 

J2 i 

2m^ r = ^d lJ U + 0{d 2 ) , (12) 
and an evolution equation for the dipole moment, 

^ = 2m g i -2F i + 0(d 2 ) . (13) 
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In both Eqs. (TT2"j) and (TT5]) we neglect terms of the order of the square of the separation 
distance d, assuming that d <C \x l \. In the RHS of (fT2l) - ffT3l) , g % and dyU are evaluated 
at the position of the center-of-mass x\ The torque produced by the forces acting on 
the two sub-particles in the RHS of (!9l)- (jT0{) is given by 



This torque will tend to align the dipole moment with (and in the same direction of) 
the gravitational field. 

The prominent feature of the equation of motion (Tl2|) is the violation of the 
equivalence principle by the dipolar particle, as a result of the fact that the particle's 
inertial mass is 2m while its gravitational mass is zero. Note that Eq. (Tl2l) has a structure 
different from the equation we obtain in the non-relativistic limit of the relativistic model 
of paper II — indeed, contrarily to the present quasi-Newtonian model, the relativistic 
model is consistent with the equivalence principle (see paper II for discussion) . While the 
inertial and passive/active gravitational masses of the dipolar particle are given here by 
Mi = 2m and M p = M a = 0, the relativistic model of paper II will have M\ = M p = 2m 
and M a = 0(c~ 2 ) in the non-relativistic limit c — > oo. This represents a fundamental 
difference between the two models. In addition, of course, since the present model is 
Newtonian it does not a priori allow (in contrast to paper II) to answer questions related 
to cosmology or the motion of relativistic particles. However, despite these differences, 
we shall recover in paper II the main characteristics of the mechanism of gravitational 
polarization. 

From Eq. ( [121) the dipolar particle is expected to accelerate slowly in a given 
gravitational field, as compared to an ordinary particle. More precisely, we find that 
the particle is not directly subject to the gravitational field, but, rather, to its gradient, 
namely the tidal gravitational field d{jU. In the potential U ~ 1/R (respectively the 
MOND analogue U ~ \nR), the acceleration is typically of the order of 1/R 3 (resp. 
1/-R 2 )@ The observational consequence is that the dark matter consisting of a fluid 
of dipole moments is necessarily cold, and even "colder" than ordinary non-relativistic 
matter. This property may be consistent with the observation of galactic structures. 
Thus, the dipolar dark matter appears as a medium whose dynamics is weakly influenced 
by the distribution of ordinary galaxies. 

On the other hand, the evolution equation ( fl3l) shows that a situation of 
equilibrium, where the distance d between the pair of particles remains constant, is 
possible. The equilibrium is realized when the internal force F l exactly compensates for 
the gravitational forcejjj] 



§ This type of coupling to the tidal gravitational field is well known; for instance it corresponds to the 
non-relativistic limit of the coupling between the spin and the Riemann curvature tensor, for particles 
with spin moving on an arbitrary background [251 126) . 

|| In Sec. [H] we shall be more precise about what is meant by equilibrium condition; see Eq. (|32[) . 




(14) 



F { = m g* + O (d 2 ) . 



(15) 
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Note that here as everywhere else, g l represents the total gravitational field, sum of the 
contributions due to the ordinary masses and the polarization masses. Because F l is 
proportional to d l , Eq. ( fiTi) . we see that when the equilibrium holds, the dipole moment 
is in the direction of the gravitational field, ir l = md % oc g l . The polarization vector 
IP = n n 1 is aligned with the gravitational field and the medium is polarized. Hence, the 
equilibrium condition ffl5|) provides a mechanism for verifying the crucial equation ([3]). 

Let us now get information on the susceptibility coefficient \ as a function of 
g = \g l \. From Eq. ([3]) and the relation U = nmd [see (JHI) — (EJ)] , we get 



The condition ffTS]) implies that d is necessarily a function of g, say d = d(g), obtained by 
inversion of the relation F(d) = mg or, rather, since as we have seen the force should 
also depend on g, of the relation F(d,g) = mg. Thus, \ is a certain function of g, 
depending on the properties of the internal force F l , and we are able, in principle, to 
relate the MOND function /i = 1 + x to the internal structure of the dipolar particles. 
We obtain 



Notice that such function is expected to be a complicated function of g, because it 
is made of the inverse of F(d,g) = mg, and especially because it depends on the 
spatial distribution of the dipole moments, characterized by their number density 
n. The distribution of n is determined by the gravitational field via the equation 
of motion f[T!?l) . together with the Eulerian continuity equation d t n + di(nv l ) = 0, 
where v l = dx l /dt. However, as we have seen, the motion of the dipolar particle is 
sensitive in first approximation only to the tidal gravitational field. Thus, a reasonable 
approximation is probably to consider that the velocity field v % remains small, hence the 
number density n is nearly constant and uniform. In the following we shall neglect the 
tidal gravitational fields, and shall treat n as a constant. 

5. Internal force law 

In the limit g — > 0, we enter the deep MOND regime — a non-linear regime characterized 
by fi = g/ao + 0(g 2 )^ Comparing with Eq. ( flTI) . we deduce that the dipole separation 
d should behave in terms of the gravitational field in the MOND regime like 



(recall that n is assumed to be constant); this means, using F = mg, that the force 



X = —A-nGmn — . 

9 



(16) 




(17) 




% For simplicity we assume a power-law expansion when g — > 0. 
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must be dominantly proportional to d. More precisely, we find that the force law F(d) 
that is necessary to account for the MOND phenomenology is 

Ait Gmn 



F(d) = \ixGm 2 nd 



1 + 



a 



■d + O (d 2 ) 



(19) 



Interestingly, this force becomes weaker when the particles constituting the dipole 
moment get closer. As we see from (I19p . the MOND acceleration scale a happens 
to parametrize, in this model, the expansion of the internal force at short distances, i.e. 
d — *■ (in the regime where g — *■ 0). 

Equation (fl9~|) represents the value of the force at equilibrium, i.e. when (fT5l) is 
satisfied, and it depends on the number n of particles. However the internal force F l 
itself — not necessarily at equilibrium — is defined by the equations of motion (jHJ)- 
( Tl0|) for a single dipole moment without reference to n, and in this sense is intrinsic 
to the dipole moment. Nevertheless it seems unusual that the equilibrium force (TT9T) 
should depend on the surrounding density n of the medium. In Sec. [6] we shall show 
that this force is actually the one corresponding to an harmonic oscillator describing the 
oscillations of a "gravitational plasma" at its natural plasma frequency. The dependence 
on n of Eq. (fl9~l) will then appear to be that involved into the usual expression [22] of 
the plasma frequency, as given by Eq. ( [291) below. 

In addition, we want to recover the usual Newtonian gravity when g ^> a , and 
we see from Eq. ( [171) that it suffices that d(g)/g — > when g — > oo. There is a large 
number of possibilities; many force laws F(d) do it in practice. For instance we can 
adopt d oc g 1 ~ e with any e > 0, which corresponds to the power force law F(d) oc d a 
with a = 1/(1 — e), and we see that any powers a, except those with < a < 1, are 
possible. In the discussion below we shall choose e > 1 in order to ensure that d —>■ in 
the Newtonian regime. The particular case e = 3/2 gives back the 1/d 2 type force law. 
In this case we can adopt for the susceptibility function 

3/2 r 



X 



OqV 

9 J 



1 + I - 

9 



This corresponds to the dipolar separation 

1/2 

d 



«0 



4:71 Gmn 



and internal force 
F = 



ma 



Oo 

9 
a 



1 + - 



0(d 2 ) 



(20) 



(21) 



(22) 



A71 G mnd' 

Notice that because the expression (1221 is positive, it represents a Coulombian force, 
i.e. attractive between unlike masses. 

Since d tends to zero in both the Newtonian and MOND regimes, we see that the 
function d — > F(d) is actually two- valued. We have already noticed that F depends 
not only on d but also on g; the expression of the force (|T9|) is valid in the MOND 
regime where g — > 0, while Eq. ( l22l) holds in the Newtonian regime where g — > oo. An 
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alternative (but rather ad hoc) choice, encompassing both types of behavior, is provided 
by the susceptibility function \ = —e~ g ^ a ° . In this case we have 

d=— | e- g/a \ (23) 

and the force law F(d) at equilibrium is obtained by substituting g = F/m in the RHS 
and looking for the two- valued inverse. 




H 

2 

Figure 1. The dipolar separation distance d as a function of g. The equilibrium 
condition F — m g is satisfied. Such graph is valid when the number density n of dipole 
moments is constant. More generally, for non constant n, the physically meaningful 
analogue of d would be the polarization H = nmd (cf. Fig. 4 in paper II). 



The generic form of the distance function d(g) is illustrated in Fig. [JJ From the 
figure we comment on the physical picture one might have in mind. In the absence of 
gravitational fields, i.e. in the absence of ordinary matter, the dipole moments do not 
exist (at least classically) since their separation d is zero. Indeed, when g = we have 
d = by Eq. ATS]) . The dipolar ether does not produce any noticeable effect. Suppose 
that some external mass, made of ordinary matter, is steadily approached. The dipole 
moments start feeling a weak gravitational field and they become active. According to 
Eq. ffTSl) and Fig. [fl the dipoles open up and get aligned in the gravitational field in 
order to maintain the equilibrium ( jl~5i) . The medium is polarized, we are in the MOND 
regime, and the gravitational field is dominated by the contribution of the polarization 
masses. Further approaching the ordinary mass, the dipolar separation d eventually 
reaches a maximal value. From (fT8l) we find, approximately, 

a ~ ^ (24) 

"max r> ri i K^J 

lo7T Gran 

which is reached for g w Oq/2. If at that point we continue to increase the gravitational 
field (putting nearer the external mass), d will begin to decrease and the dipole moments 
will close up. Finally, for strong gravitational fields, the dipole moments become inactive 
again (indeed d = when g — > oo). The gravitational field is dominated by the 
contribution of the ordinary matter, and we recover the Newtonian regime. 



6. Gravitational plasma 

The medium of dipolar dark matter described above can also be interpreted, in an 
alternative formulation, as a polarizable "gravitational plasma" consisting of the two 
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species of particles (mi,m g ) = (m, ±m). Such formulation will essentially be equivalent 
to the previous one, but is simpler and more appealing physically. 

Suppose that the two types of particles (m, ±m) are in equal numbers so that the 
plasma is globally neutral. There is an equilibrium point where the plasma is locally 
neutral, so the number densities of the two particle species are equal, at each point, to 
some common constant and uniform value n. As it stands the equilibrium is unstable 
because the gravitational force between unlike masses is repulsive. To ensure a stable 
equilibrium we must postulate like in Eqs. (|9l)-(fT0l) some restoring non-gravitational 
force F l , acting between the masses m g = ±m, and superseding the gravitational force. 
We introduce some associated internal field f l such that 

r = -mf\ (25) 

and assert that this field obeys a Gauss law in the non-relativistic limit, 

ft/ < = -^^(n + -n_). (26) 

X 

The number densities of the particles species (m, ±m) are denoted by n±. We have 
introduced here the susceptibility coefficient \ t° represent a dimensionless coupling 
constant characterizing the internal interaction. Notice that since x is negative, Eq. (jSJ), 
the force law ff26l) is attractive between unlike particles and repulsive between like 
ones. Furthermore, supposing that the plasma is bathed by an external gravitational 
field g l } constant and uniform over some region of consideration^] we expect that the 
coupling constant should reflect the presence of this gravitational field, and we assume 
a dependence on its norm g, namely \ = x{q)- 

Let x\ and x % _ (in short x±) be the displacement vectors of the masses from the 
equilibrium position characterized by the density n. The particles are accelerated by 
the internal field f l and by the applied external gravitational field g l . The equations of 
motion have already been given in Eqs. (l9l)- (|T0l) and read now 



J2 i 

m-^ = ±m(f + g t). (27) 

Consider a small departure from equilibrium, corresponding to small displacement 
vectors x l ±. The density perturbations are given by n± = n (1 — diX l ±) to first order 
in the displacements x l ±. Using Eq. (f26|) we readily integrate for f % and inject the 
solution into ( J27j) . In this way we find that x\ + x l _ = (the centre of inertial masses 
is at rest — neglecting tidal fields), together with the following harmonic oscillator for 
n l = m d l where d % = x\ — x l _ , 



dt 2 



+ u 2 V = 2m g i . (28) 



+ We adopt the frame associated with the plasma's equilibrium configuration. The gravitational field 
g l is defined in that frame. 
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Actually this computation is the classic derivation of the plasma frequency [22J which 
is found for the case at hands to be 



, 87r G m n , „ 

c i / — . (29) 

The frequency depends on the density n of the plasma at equilibrium and on the strength 
of the internal interaction which is encoded into the coupling constant \- The solution 
we get for the internal field is 

2 2 

/* = -— it* F' = — tt\ (30) 

J 2m 2 ' V ; 

so we can check that Eq. (128]) is equivalent to our previous equation ({TBI) (recall that 

here g i is considered to be constant and uniform). The particles oscillate around some 

non-zero mean position that is determined by the ambiant gravitational field as 

<**> = ^. (31) 

The mean value of the force is 

(F i )=mg i , (32) 

and coincides with the equilibrium condition postulated in Eq. ( fl5l) . 

This result ( I3T1) is classical but can easily be recovered from standard quantization 
of the harmonic oscillator (1281) . The spectrum of "plasmons" is discrete (k G N) with 
energy levels shifted by the gravitational field |27J, 

„ A 3\ mq 2 , . 

E>=[k +i )hu,-JL. (33) 

The eigenstates are \ipp lP2P3 ) = \4>pi) \4>p 2 ) I ^3) where k = P1+P2+P3', the eigenvalues have 
degeneracy (k + l)(k + 2)/2. The one- dimensional eigenstate function Pi (vr l ) = (7r*|0 p .) 
is of the form ^> Vi [^ 1 ~ ^t9 l ) where (H p being the Hermite polynomial) 

Vp(p) = ^ f^ L ) 1/4 ^(v / I/ , ) e '^ 2 - (34) 
22^! \7rnm/ v am 

The expectation value (tt 1 ) = (V , pip 2 p3l 7r? lV'pip2P3) does not vanish due to the presence of 
the gravitational field; it can be computed using f_™dinr \(p p {^ — ^g)\ 2 = ^jrg and 
the result agrees with ( 13TI) . 

We analyze the effect of the mean dipole moment vector (!3~T!) on the equation for 
the gravitational field g % (supposed now to be generated by external sources). More 
precisely we argue that the mean value is really a quantum expectation value as we 
have just proved, so that by adopting a "semi-classical" approach this expectation value 
(tt 1 ) should have a direct effect as the source for the classical gravitational field g l . The 
polarization readily follows from the expression of the plasma frequency (|29|) as 

W = n(^) = -^g\ (35) 

so that we recover exactly Eq. ([3]). Note that the constant density n of the plasma at 
equilibrium is equal to the density of dipole moments. The polarization is automatically 
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proportional to the gravitational field. This is an attractive feature of the present 
formulation based on a gravitational plasma; there is no need to invoke a mechanism to 
align IP with g l . Following the same reasoning as in Sec. EJ the polarization ( [351) gives 
rize to the density of polarization masses (jl]), which is added to the RHS of the Poisson 
law, Eq. (jHJ). The MOND equation (JTJ readily follows. 

Let us comment on the field equation fl26|) for the supposedly fundamental internal 
interaction f % . In the absence of gravity we have \ = — 1, hence Eq. fl26|) becomes 
that of a negative Poisson equation corresponding to a negative Newtonian force (with 
"anti-gravitational" constant —G). On the other hand, in strong gravity (g/a — > oo) 
we have x ~ * an d we see that the strength of the internal field f % becomes infinite. 
Therefore the plasma gets locked in its undisturbed equilibrium state for which 7r* = 
(strictly speaking this is true if (7X — ^ when g — > oo). In this limit the dipolar medium 
is inactive; there is no induced polarization and the Newtonian law holds. 

Note that there is clearly some amount of fine tuning in the present model (and 
the one of paper II). Namely the function xio) is n °t determined from first principles 
within the model but is tuned to agree with the phenomenology of MOND. In particular 
%(0) = —1 is not "explained" in the model but comes from astronomical observations. 
Similarly = is imposed in order to recover the Newtonian regime where it is 

observed to be valid. 

To conclude, the phenomenology of MOND suggests the existence of a polarization 
mechanism at work at the scale of galactic structures, and which could be viewed as 
the gravitational analogue of the electric polarization of a dielectric material. The 
dark matter would consist of the polarization masses associated with some gravitational 
dipole moments aligned with the gravitational field of ordinary masses. We find that 
the properties of the dipolar dark matter are governed by the internal non-gravitational 
force linking together the constituents of the dipolar medium. The formulation of this 
medium in terms of an oscillating gravitational plasma polarized in the gravitational 
field is particular attractive. In paper II we shall show how the notion of dipolar particle 
can be made compatible with the framework of general relativity. 
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